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Abstract 

A number of numerical simulations of lattice gauge theory have indicated a low 
mass of strange quark in 100 MeV range at the scale of = 2 GeV. In the 
unquenched case, which is improved over the simulation in the quenched approxi- 
mation by the inclusion of u and d sea quark effects, one sees a further downward 
trend. Here the fermion mass spectrum of the Fritzsch texture is recalculated. 
In a single step supersymmetric GUT with Mx ~ 10^^ GeV such values of the 
strange quark mass can be obtained for low values of tan /?. Experimental numbers 
^poie _ -^Y^ -|- g QeV and 4.1 < mi,{mb) < 4.4 GeV are used in this study. Since 
the scenario is supersymmetric, gaugino loop diagrams contribute to the masses 
in addition to usual tree level Yukawa contributions. Upper bound of the mixing 
parameter Vet is taken at 0.045 
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The strange quark mass has been traditionally calculated using the cur- 
rent algebra mass ratio 

=12.6 ±0.5 (1) 

rriu + TUd 

Equation (|l]) is evaluated using values of (m^ + m^) which are reported in 
calculations of QCD finite energy sum-rules (FESR). At the two-loop level 
of perturbative QCD calculations which include non-perturbative corrections 
up to dimension six one has the result 

(m„ + mrf)( 1 GeV) = 15.5 ±2.0 MeV (2) 

For as{mz) = 0.118 results of (|l]) and (0) together lead to 

m,(l GeV) = 195 ± 28 MeV or m,(2 GeV) = 150 ± 21 MeV. (3) 

The ratio 

m,(2 GeV)/m,(l GeV) = 0.769 for ^.(m^) = 0.118 (4) 

can be obtained by solving renormalization group equations^]. A systematic 
uncertainty in this result remains in reconstruction of so called 'spectral 
function' from experimental data of resonances. When a different functional 
form of the resonance is adopted, and three loop order perturbative QCD 
theory is used one obtains 0] 

(m„ -t- ma) (1 GeV) = 12.0 ± 2.5 MeV (5) 

With dl) and (|ip one gets 

m,(lGeV) = 151 ± 32 MeV or m,(2GeV) = 116 ± 24 MeV. (6) 

Again this translation from the scale of 1 GeV to the scale of 2 GeV is 
obtained for the case ag = 0.118. It has been remarked in Ref. [^] that 
it is indeed difficult to account for vacuum fluctuations, or sea quark effects 
generated by quarks of small masses in perturbative QCD calculations. Thus, 
numerical simulations of strange quark mass on a lattice becomes rather 
attractive, especially if the simulation includes virtual light quark loop effects. 
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Up and down type quarks differ only in the U{l)em quantum numbers 
in an effective theory where the gauge symmetry is SU{3)c x U{l)em- Lat- 
tice calculations in current literature have neglected effects of f/(l)em which 
distinguishes up quarks from down quarks. Let us note that we are describ- 
ing the lattice in terms of a theory at the scale of a few GeVs where light 
quark masses are to be described in terms of observables relevant to their 
own scales, which are meson masses and decay constants. Thus, lattice sim- 
ulation determines m^, '^"'^"^^ and the lattice spacing a using three hadronic 
observables. They can be chosen, for example, Mj^, M^*, and Due to the 
structure of equations which are to be fitted, the scale a can also be taken as 
a function of some other observable, for example, it may be chosen as a{Mn) 
or a(MA) etc. The result depend on the choice of the observable that fits 
the lattice spacing. The best choice would be the one which has minimum 
experimental uncertainty and the best result would be a clever weighted av- 
erage of results from various choices. A test of the simulation is obviously to 
see whether results from various observables are statistically consistent with 
each other. 

Next question is how do we describe quark masses when the theory is 
living on a discritized lattice. Various definitions or formalisms of quark 
masses on a lattice have been suggested. Ref.p uses the definition in terms 
of hopping parameter k of the lattice 

a rribare = log (1 + (1/2/t - 1/2/tc)), (7) 

for a Wilson-like fermion. In the continuum limit we have a ^ 0, and there 
one gets the hopping parameter k = Kc = 1/8. A smaller hopping parameter 
makes the lattice more sticky, and fermions remain on lattice points for a 
longer time. This make them look as if they were more massive. 

There are various other formalisms of defining the mass of fermions on the 
lattice such as stagerred fermions or domain wall fermions. Most calculations, 
however, use the Wilson action for various definitions of the fermion mass. 
Next, to compare the result with experiment, one has to calculate the MS 
mass at a scale /i starting from the lattice estimate of the bare mass (|^) 
using, for example, the mass renormalization constant Zm(/^) relating the 
lattice regularization scheme to the continuum regularization scheme. The 
lattice regularization prescription is given in Ref[0. The Zm constant for 
various formalisms such as Wilson-like or Staggered are given in table^ of 
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[^]. Final results of the physical quark mass for various definitions of the 
fermion on a lattice differ 0{a) among each other and one expects to get the 
same result of the physical quark mass in the continuum limit when a — > 0. 

Beyond the minimal lattice simulation of light quark masses using the 
heavy quark effective theory, the next step would be to incorporate sea quark 
effects. From the conservation of energy it can be understood that it is easiest 
to produce hghtest of the quarks virtually. Indeed such simulations have 
been performed. They are termed nf = 2 unquenched lattice simulations. 
The detailed processes of numerical simulations are described by respective 
authors. However, we have summarized the results of recent studies are in 
table I in the order: (A)i, (B)i, (C)!!, (D)|lll, (E)i, (F)|l, and 
(0)0. 



reference 


quenched 


dynamical 


(1/ a)callihration 


nis{MeV) 


A 


yes 




nip 


143 ±6 & 115 ±2 


B 


yes 




niK* 


130 ± 20 


C 


yes 




niK* 


122 ± 20 


D 


yes 




niK* 


111 ± 12 


E 


yes 




nip 


110 ±31 


F 


yes 




nip 


108 ±4 


G 


yes 




IP — IS splitting 


95 ±16 


A 




yes 


nip 


70 & 80 


E 




yes 


nip 


68 ±19 


G 




yes 


IP — IS splitting 


54-92 



Table 1: Ref. G uses IP-IS splitting of the charmonium system to calibrate 
(1/a). Reference A quotes two different results for two sets of parameters. 
All results are at the scale fi = 2 GeV 

On the experimental front the bottom quark mass is in the range 

4.1 < mb{mb) < AA GeV (8) 
according to the review of particle physics (PDG) tables[|T5. Theoretically, 



one re-expresses the bottom quark mass in terms of parameters of the Min- 
imal Supersymmetric Standard Model(MSSM). The tree level contribution 
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which is related straight to the Yukawa texture, and the one loop contribu- 
tion due to the dominant gaugino loop can be accounted individually. Then 
one can write down the relation 



rrib = ml''''''''' + mf^'^^ 

Vf 82 tan/5 TTig /i 

= h,-j^cosP + m,-g,^^-^. (9) 

Here rrig is the gluino mass /i is the fi parameter and ^e// is averaged su- 
per symmetry breaking mass scale. This paper discusses a scenario where 
the first term of the RHS of (|^) comes from diaginalizing a Fritzsch Yukawa 
texture. The second term can be estimated to be around ±2 GeV. In the su- 
persymmetric case it will be satisfactory if the Fritzsch Yukawa contribution 
is in the range 

2.1 < m*"^*""" < 6.4 GeV (10) 

Next question is concerning tan /3. Supersymmetry, together with the gauge 
quantum number structure of fermions demands that at least two Higgs dou- 
blets are necessary. Hence the ratio of the VEVs of two Higgs doublets is an 
unavoidable parameter given the value of the effective four-Fermi coupling 
Vp. There are perturbative bounds on tan j3 in the context of grand unified 
theories (They can be extended to supersymmetric theories without grand 
unification if MSSM is valid up to a certain high scale, say 10^^ GeV). In 
practice the very low- valued regions of the parameter space for tan (3 are for- 



bidden from perturbative considerations. See |T6[ for example. Furthermore 
high values of tan/3 have constraints from charge and color breaking[17, |l^] . 
This is especially true if Yukawa couplings are at the fixed point region. The 
intermediate regions of tan (3 are definitely allowed. To make a safe case let 
us choose the range for the purpose of this paper 

tan/? = 2-30. (11) 

Now let us focus on the texture. It has been noted that the quark mixing 
angle Vus, which is a dimension-less quantity, can be thought of as a ratio of 
the mass scales of flavor symmetry breaking. These symmetries lead to the 
mass hierarchy between families. Phenomenologically of course, the ratio of 
the masses of the first and the second generation satisfies well the relation 



t.ne.-0 (12) 
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If there are two Higgs doublets instead, ( [I^ ) remains untouched as the ratio 
of the VEVs of the doublets cancel in the ratio on the RHS. Thus it cannot 
feel tan /3. 

Suppose in a two generation case rotation angles of the up and the down 
sectors are 9u and 9d- Then the combined quark mixing matrix V = OjD^ 
will give 6c = Ou i: 6d- This observation plays a role in the Fritzsch mass 
matrices. Fritzsch mass matrices can be thought of as a set of mass matrices 
which generalizes (|12D to the following form[|19|, ^ 



9c = 9d± 9u 



tan 




(13) 



Where rrii are eigenvalues of Fritzsch mass matrices. In the three generation 
case Fritzsch textures for up and down sectors are given by 

/ ae*^ \ / Ae'^ 

Md = 



M, 



u 



ae 











(14) 



The phases r, R, r', R', h, H, h', H', q, Q can be absorbed in the redefinition 
of quark fields and individual mass matrices can be made real. However 
the weak charge changing current consists of (couples to) gauge eigenstates. 
Furthermore the structure of the weak charge changing current is ijj^'^^ip^. 
Consequently the weak mixing matrix must contain some combination of 
phases. It can be shown that residual phases of the weak mixing matrix are 



H) 



a = (r-R)~{h 
r = {r-R)-{h' 

when the weak mixing matrix is expressed as 

/I 



ih' - H') + {q-Q) 



O 



u 





Vo 



o 



D 



(15) 



(16) 







In (|T6D Ou and Od diagonalizes Aiu and M.d which are precisely those in 
(|np but in the limit when all the phases vanish. Which are very simply 



(17) 





/O 


a 






/O 


A 




Mu = 


a 





;) 


; Md = 


A 





i) 




lo 


h 






lo 


B 





6 



Let us consider the limit tUu << rric « rrit and << << rrib. We can 
approximately re-write (^) as following 



Ml 




M 



D 




(18) 



Now we easily see that M.u can be diagonalized by the rotation 



Of/ = I -sin^]^ cosO'l 
1 



1 

cos sin 

-sin^^ cos^^ 



(19) 



where tan 9'^ = and tan O2 = -y/^- M. d can be diagonalized similarly. 
In the limit of strongly hierarchical eigenvalues one may use the approxima- 
tion 



cos^r 



COS0° 



Furthermore let us define the quantities 



V rrir. 




1^1 




1^2 




(20) 



(21) 



Then the mixing matrix (|16|) takes the following form 



/ 1 

-^i + Vie'" 



/ii(z/2e'" - yU2e' 



/i2e' 



(22) 



A detailed derivation of these relations are given in Ref|^. A cancelation 
among two terms in the expression for Vcb in (^2D is needed. To achieve this 
one can choose 



TT 



cr ~ 7- ~ — — this gives Vcb ~ W W — • 



(23) 



It is easy to check that (|23|) makes the top quark mass too light to be exper- 
imentally true. Thus, it is worth asking the question whether if the Fritzsch 
texture were valid at the GUT scale instead, in other words, if the flavor 
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symmetries were exact only above the GUT scale, could a miracle of renor- 
malization group evolution of the masses and mixing angles make the Fritzsch 
relations valid at low energyP^. Here we study their idea. 

The renormalization of the full 3x3 complex Yukawa matrices and their 
renormalization up to the GUT scale Mx = 10^^ GeV is what follows. Let us 
set our notations of mixing angles, the non-removable phase and eigenvalues 
of the Yukawa matrices. We adopt the parameterization 



S1C2C3 - ciS2e*'^ C1C2C3 + siS2e"^ C2S3 

V ~'SiS3 — C1S3 C3 



(24) 



There is a detailed proof in Ref . ||23| that in this parameterization eigenvalues 
Ui of the Yukawa textures, three CKM mixing angles and the CP violating 
phase (f) satisfy the renormalization group equations 



dt 


= , 








9 d 

167rV 
dt 


In tan 9i 


= -yl sin^ ^3 , 






9 d 

levrV 
dt 


In tan 6*2 


= -yl sin^ 6*3 , 






9 d 

levrV 
dt 


In tan 9^ 


2 2 






9 d 
dt 


In l/„ = - 


-dtgl + 3y,2 + yl cos' ^2 


sin' 03 




9 d 
dt 


In ?/c = - 


-cl'(7,2 + 3y,2 + sin^ ^2 


sin' 63 




d 
levrV 
dt 


In l/i = - 


-dtg} + 6i/2 + y^cos^^s 






9 d 
dt 


In = - 








^a 2 d 
Wit — 

dt 


In 1/s = - 


-cfgj + cos' 61 sin' ^3 


+ 3z/b' 




r, d 

167rV 
dt 


In ?/b = - 


-cf (7,2 + cos' 03 + 6?/,' 


+ 2/^ 




9 d 
dt 


In 1/e = - 


<9l + ^yl + yl , 
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levr^^ In = -0^(72 ^3y2 + 4y2 . (25) 



Multipliers c*- of gauge couplings in (pTsP are well known. They can be found 
in ^3|]. We have solved these one-loop equations numerically using Math- 
ematica NDSolve subroutine. The flow chart follows this line. Taking all 
experimentally possible values of the eigenvalues but only central values of 
the angles at low energy, we have evolved the set to the GUT scale using 
At the GUT scale predictions for Vcb Vus and Vub are calculated as- 



suming that Fritzsch relations are valid only at the GUT scale and beyond. 
While translating predictions of CKM entries back to low energy using ([T8|), 
we have used exact values of the angles not central values. Thereafter we 
have checked whether each individual value of masses and mixings remain 
within experimentally allowed ranges. For the strange quark mass values 
quoted in table|l| are used. For all other masses and mixings the experimen- 
tal values are taken from the review of particle physics Our results are 
given in tableEl 



Ots 


tan (5 


m,(2 GeV) 


0.118 


2 


59.90 MeV 


0.118 


10 


61.52 MeV 


0.118 


20 


63.05 MeV 


0.118 


30 


66.90 MeV 



Table 2: Our results are quoted for mf = 173 GeV. All other masses and 
mixings remain within the ranges quoted by the Review of Particle Physics. 

In conclusion, implications of results of Uf = 2 unquenched lattice simu- 
lations of the strange quark mass in the context of the Fritzsch texture are 
studied. Previous calculations in this line exist in the literature. We have 
two new aspects. Because the combined effect of charge and color breaking 
and perturbative unitarity of the Yukawa couplings may rule out the large 
tan/5 scenario we have studied the low tan/3 scenario. Moreover, we have 
included corrections of supersymmetric origin in the study. 
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Supersymmetric corrections to the bottom quark mass and the low tan f3 
scenario goes hand in hand. This is in the sense that in the low tan f] regime 
Fritzsch texture demands a large Yukawa contribution to the bottom quark 
mass. This is partially canceled by the supersymmetric loop corrections. 
Thus, the original Fritzsch texture is consistent with experimental data if it 
holds at the GUT scale. The strange quark mass emerges in the range 60-70 
MeV at /i = 2 GeV for the central values of = 0.118 and mf"'*^ = 173 GeV. 
This range is consistent with nj = 2 sea quark effect improved (unquenched) 
lattice simulations of the strange quark mass. 

We thank S. Gottlieb for discussions and M. S. Berger for lending his 
Fortran routine to calculate the rji scaling factors of Ref.0]. This research 
was supported by U.S Department of Energy under the grant number DE- 
FG02-91ER40661. 
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